BOUNDARY QUASI-ORTHOGONALITY AND SHARP INCLUSION 
BOUNDS FOR LARGE DIRICHLET EIGENVALUES 



A. H. BARNETT AND ANDREW HASSELL 

Abstract. We study eigenfunctions 4>j and eigenvalues Ej of the Dirichlet Laplacian on a 
bounded domain Q C K" with piecewise smooth boundary. We bound the distance between an 
arbitrary parameter E > and the spectrum {Ej} in terms of the boundary L^-norm of a normalized 
trial solution u of the Helmholtz equation (A + E)u = 0. We also bound the L^-norm of the error 
of this trial solution from an eigenf unction. Both of these results are sharp up to constants, hold 
for all E greater than a small constant, and improve upon the best-known bounds of Moler— Payne 
by a factor of the wavenumber \/E. One application is to the solution of eigenvalue problems at 
high frequency, via, for example, the method of particular solutions. In the case of planar, strictly 
star-shaped domains we give an inclusion bound where the constant is also sharp. We give explicit 
constants in the theorems, and show a numerical example where an eigenvalue around the 2500th 
is computed to 14 digits of relative accuracy. The proof makes use of a new quasi-orthogonality 
property of the boundary normal derivatives of the eigenmodes (Theorem 11.31 below) , of interest in 
its own right. Namely, the operator norm of the sum of rank 1 operators dn<j)j {dn(f>j , ■) over all Ej 
in a spectral window of width V~E — a sum with about terms — is at most a constant 

factor (independent of E) larger than the operator norm of any one individual term. 

1. Introduction and main results. The computation of eigenvalues and eigen- 
modes of Euclidean domains is a classical problem (in two dimensions this is the 
'drum problem', reviewed in [211134) 1 with a wealth of applications to engineering and 
physics, including acoustic, electromagnetic and optical cavity and resonator design, 
micro-lasers , and data analysis [30] . It also has continued interest in mathemati- 
cal community in the areas of quantum chaos |37[ |3] and spectral geometry |18| . Let 
be a sequence of orthonormal eigenfunctions and Ej the respective eigenvalues 
(0 < El < E2 < E3 < ■ ■ ■ counting muhiplicities) of -A, where A := Em=i d'^/dx^, 
is the Laplacian in a bounded domain Q G M", n > 2, with Dirichlet boundary 
condition. That is, (jjj satisfies 

{A + Ej)(j>j^Q inn (1.1) 

(t)j = on dn (1.2) 

||0j||L2(f2) = 1. (1.3) 

We will call the spectrum a :— {Ej}°°^i. Many of the applications mentioned demand 
high frequencies, that is, mode numbers j from 10^ to as high as 10^. Efficient 
solution of the problem thus requires specialized numerical approaches that scale 
with wavenumber better than conventional discretization methods. 

The goal of this paper is to bound the errors of approximate eigenvalues and 
eigenfunctions computed using trial functions that satisfy exactly the homogeneous 
Helmholtz equation in il. As we will review below, such computational methods have 
proven very powerful. Recently one of the authors [4] improved upon the classical 
eigenvalue bound of Moler-Payne |24j by a factor of the wavenumber; however, this 
result has limited utility since it applies only to Helmholtz parameters lying in neigh- 
borhoods of a of unknown size. In the present paper we go well beyond this result 
by giving new theorems, which i) hold for all Helmholtz parameters (greater than an 
0(1) constant), ii) retain the improved high-frequency asymptotic behavior of [41 and 
show that this behavior is sharp, and iii) improve upon the best-known eigenfunction 
estimates, again by a factor of the wavenumber. To achieve this we make use of a 
new form of quasi-orthogonality of the eigenfunctions on the boundary, Theorem II. 31 
of independent interest. 
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Fig. 1.1. Tension t[umin] versus energy E for the domain shown on the right. Umin is the optimal 
trial Helmholtz solution lying in the span of a numerical basis set (see Section^, a) Low freguency, 
showing the minima corresponding to the lowest 20 Dirichlet eigenvalues, b) Medium-high frequency, 
showing a similar interval starting at eigenvalue number j 2552/ note the new vertical scale, c) 
Density plot of eigenfunction rf>j x itmin corresponding to the eigenvalue Ej = 10005.02135797 ■ ■ • 
shown by the dot in b) (black indicates large values of \4>j\'^ , white zero). 



Before presenting our results, we need to review some known inclusion bounds 
and their importance for applications. Given an energy parametei0 E > 0, let u be a 
non-trivial solution to the homogeneous Helmholtz equation {A + E)u — in fl with 
no imposed boundary condition, and define its boundary error norm (or 'tension') 

\W\\L^n) 

Clearly, t[u] =0 implies that E is an eigenvalue. It is reasonable to expect that if t[v] 
is small for some Helmholtz solution m, then E is close to an eigenvalue. Moler-Payne 
[24] (building upon [16]) quantified this: there is a constant Cmp depending only on 
the domain, such that 

d{E,a) < C^pEt[u] , (1.5) 

where d{E,a) :— minj \Ej — E\ denotes the distance of E from the spectrum. 

An important application is to solving (jl.ip - (jl.3l) via global approximation meth- 
ods, including the method of particular solutions (MPS) [9j IJ. One writes a trial 
eigenmode u = X]n=i Cn^n via basis functions which are closed-form Helmholtz 
solutions in D, but which need not satisfy any particular boundary condition. By ad- 
justing the coefficients c„ (via a generalized eigenvalue :4! or singular value problem 
[5j) one may minimize t[u\ at fixed E; by repeating this in a search for E values where 
the minimum t[u\ is very small, as illustrated by Fig. Il.lb and b, one may then locate 
approximate eigenvalues whose error is bounded above by (jl.5|) . (This is sometimes 
called the method of a priori-a posteriori inequalities [2h Sec. 16].) 

Due to the work of Betcke-Trefethen [9] and others, such methods have enjoyed a 
recent revival, at least in n = 2, due to their high (often spectral) accuracy and their 
efficiency at high frequency when compared to direct discretization methods such as 
finite elements. For example, in various domains, 14 digits may be achieved in double 



^The Helmholtz parameter E may be interpreted as energy, or as the square of frequency, de- 
pending on the application. 
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precision arithmetic [5] , and with an MPS variant known as the scahng method, tens 
of thousands of eigenmodes as high as j ~ 10^ have been computed [3J [S]. (There 
are also successful variants [UJ |TS] by Descloux-Tolley, Driscoll, and others, in which 
subdomains are used, which we will not pursue here.) 

If we instead interpret u as the solution error for an interior Helmholtz boundary- 
value problem (solved, for instance, via MPS or boundary integral methods), then 
(|1.5|) states that the interior error is controlled by the boundary error; this aids the 
numerical analysis of such problems [211 [B]. Similar estimates (which, however, rely on 
impedance boundary conditions) enable the analysis of least-squares non-polynomial 
finite element methods [25, Thm 3.1]. Improving such estimates could thus be of 
general benefit for the numerical solution of Helmholtz problems. 

Recently one of the authors [4, observed numerical evidence that (|1.5I) is not sharp 
for large E, and showed that there is a constant Cb depending only on {I, such that, 
for each e > 0, 

d{E,a) < C^{l + e)VEt[u] (1.6) 

holds whenever E lies in some open, possibly disconnected, subset of the real axis 
containing a. This is an improvement over ()1.5|) by a factor of the wavenumber v^, 
which in problems of interest can be as high as 10^. However, since the proof relied 
on analytic perturbation in the parameter E, there was no knowledge about the size 
of this (e-dependent) subset, hence no way to know in a given practical situation 
whether the error bound holds. The point of the present work is then to remedy this 
problem by removing any restriction to an unknown subset, and also to extend the 
^/E improvement to bounds on approximate eigenfunctions. 

We assume the domain fl C R" has unit area (or volume for n > 2), and obeys 
the following rather weak geometric condition. 

Condition 1. The domain C K" is bounded, with piecewise smooth boundary 
in the sense of Zelditch-Zworski 1381. This means that 51 is given by an intersection 

N 

0= f|{x| ,A(x) >0}, 

i=l 

where the fi are smooth functions defined on a neighborhood of CI such that 

• V/, 7^ on the set {fi = 0}, 

• {fi — fj ^ 0} '"^ embedded submanifold ofR", 1 < i < j < N, and 

• fl is locally Lipschitz, i.e. for any boundary point xq G dQ, there is a Eu- 
clidean coordinate system zi , . . . , z„ and a Lipschitz function k of n — \ vari- 
ables such that in some neighborhood ofn.^, we have 

(917 = {z„ = fc(zi,...,z„_i)}. (1.7) 

Our main result on eigenvalue inclusion is the following. 

Theorem 1.1. Let fl C M" be a domain satisfying condition [TJ Then there 
are constants C, c depending only on 51, such that the following holds. Let E > 1 
and suppose u is a non-trivial solution of (A + E)u ~ in C°°(il), with t[u] := 
l|w||L2(sn)/||"||L2(n)- Then, 



d{E,a) < C^/Et[u] , 
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(1.8) 



and for the normalized Helmholtz solution u^^^ minimizing t[u\ at the given E, 

c%/!Et[u„,„] < d{E,a) < C\/^t[u„„]. (1.9) 

Remark 1.1. The estimate states that (jl.Sp is sharp, i.e., using t[u\ alone 

one cannot localize the spectrum any more tightly than this, apart from optimizing the 
constants c and C . 

Remark 1.2. The existence of a minimizer for t[u\ follows from Lemma 
in the case that E is not a Dirichlet eigenvalue ( and is trivial when E is a Dirichlet 
eigenvalue). The lower bound on the distance to the spectrum in (jl.9p is of use when 
the numerical scheme is known to produce a good approximation to u„in. 

We will also prove the following corresponding bound on the error of the trial 
eigenfunction u, which improves by a factor \/E the previous best known result 
(Moler-Payne [21 Thm. 2]). 

Theorem 1.2. LetQ be as in Theorem \l.l\ Then there is a constant C depending 
only on VL, such that the following holds. Let E > 1, let Ej be the eigenvalue nearest 
to E, and let Ek the next nearest distinct eigenvalue. Suppose u is a solution of 
{A + E)u — in C°°{fl) with ||u||i2(-Q) = 1, and let iij be the projection of u onto the 
Ej eigenspace. Then, 

\\u'U,\\LHn) < ■ (1-10) 



Remark 1.3. The left-hand .side above is equal to sinO, where 9 is the subspace 
angle between u and the Ej eigenspace (this viewpoint is elaborated in '9, Sec. 6]). For 
example, when Ej is a simple eigenvalue, we may write \\u — 4'j\\L'^(n) — 2sin(0/2). 

Remark 1.4. This result is also sharp, in a certain sense: see Remark\J^.l\ 



To conclude the introduction, we present some key ingredients of the proofs. 
Define the boundary functions of the eigenmodes by 

^^i^s) := dn(i3j{s) sedn (1.11) 

where 9„ = n • V is the usual normal derivative. Our main tools will be two theorems 
stating that boundary functions ipj lying close in eigenvalue are almost orthogonal. 
The first is the following new result which we prove in Section [51 

Theorem 1.3 (spectral window quasi-orthogonality) . Let C K" be a domain 
satisfying Condition\^ There exists a constant Co depending only on Q. such that the 
operator norm bound 

V < CnE (1.12) 

\Ej~E\<E^/-^ 

holds for all E > 1. (Here, (•, •) denotes the inner product in L^(dfl).) 

Remark 1.5. By Weyl's Law Ch. 11] there are 0{E'^^^^^/^) terms in the 
above .sum. Since each term already has norm > cE \28\ [131, theorem expresses 
essentially complete mutual orthogonality, up to a constant. Only the scaling of the 
window width with E is important: the theorem also holds for a window \Ej — E\ < 
cE^^^ for any fixed c (Cn will then depend on c as well as fl). On the other hand, one 
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could not expect it to hold over a spectral window of width 0{E^) for /3 > 1/2, since 
the boundary functions are approximately band-limited to spatial wavenumber E^/'^ 
and thus no more than 0(ii'("~^)/^) of them could be orthogonal on the boundary. 

The second result is a pairwise estimate on the inner product of boundary func- 
tions lying close in eigenvalue, with respect to a special inner product: (Here, x(s) 
refers to the location of boundary point s relative to a fixed origin, which may or may 
not be inside fi.) 

Theorem 1.4 (pairwise quasi-orthogonality) . LetQ, C M" be a bounded Lipschitz 
domain, and let S := ^ sup^^gj^ ||x|| . Then, for all i,j > 1, 



(x(s)-n(s))V,(s)^j(s)ds-2£;,% < S'{E,-EjY (1.13) 

9n 

Remark 1.6. This theorem was proved by the first-named author in Ap- 
pendix B]. It may be viewed as an off-diagonal generalization of a theorem of Rellich 
\28^ which gives the i = j case. The boundary weight x-n ( also known as the Morawetz 
multiplier) is the only one known that gives quadratic growth away the diagonal yet 
also gives non-zero diagonal elements. 

Note that neither of the above quasi-orthogonality theorems implies the other. 
We also note that Backer et al. derived a completeness property of the boundary 
functions in a (smoothed) spectral window [2j Eq. (53)], that is closely related to 
Theorem O 

After proving Theorem 11.31 we combine it with a boundary operator defined in 
Section[3]to prove the main theorems, in Section|4l In Section[5]we state and prove a 
variant of Theorem 1 1.1 1 for strictly star-shaped planar domains, which has an optimal 
constant C. This builds on Theorem 11.41 combined with the Cotlar-Stein lemma (see 
Lemma [5T^ . In the main Theorems 11.11 FOl and lOI the domain-dependent constants 
are not explicit; we discuss their explicit values in Section [HI We present a high- 
accuracy numerical example using the MPS, and sketch some of the implementation 
aspects, in Section [T] Finally, we conclude in Section [H 



2. Quasi-orthogonality in an eigenvalue window. Here we prove Theo- 
rem ll.3l using a "TT* argument". We need the fact that the upper bound ||'0j"|li2(gf2) < 
CEj on eigenmode normal derivatives, proved for example in |19', generalizes to quasi- 
modes living in an 0{E^^^) spectral window. The proof is almost the same as in [T5] . 

Lemma 2.1. Let C M" satisfy Condition[l[ Let E > 1, and let 

<!>■= E ^^■'^J (2.1) 

\Ej-E\<E^/^ 

with real coefficients cj, and — ||0|||2(52) = 1- Then, 

Wdn^Uan) < CnE (2.2) 

where the constant Co depends only on Q. 

Proof. To prove this we need the following lemma, proved in Appendix |X1 stating 
that for any piecewise smooth domain (in the sense of Condition [1]) there is a smooth 
vector field that is outgoing at each boundary point. 

Lemma 2.2. Let Q satisfy Condition]^ Then there exists a smooth vector field 
a, defined on a neighborhood ofQ, such that 



a • n > 1 
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(2.3) 



on 



almost everywhere on dQ. 

The main tool for proving Lemma 12.11 is the identity 

iD(b)d„^ = - / (/)[A, D]q} + f iD^)iA + E)(j)- f 4>D{A + E)4> (2.4) 
Jn Jn Jn 

for any first order differential operator D, which follows fro Green's 2nd identity, 
the definition of the commutator, and = 0. Choosing := a • V, where a is as 
in Lemma [2?2l we notice that the left-hand side of (|2.4p bounds the left-hand side of 
(1121), since 

/ V'l < / (a-n)V'l (2.5) 

by Condition [TJ We may now bound each of the terms on the right-hand side of (|2.4p . 
Defining Ca — sup^gn 1^(^)1: have 

mWh^n^ < Cl I ||V0||2 = --d I ^A^ = Cl J2 < (2.6) 

Jn Jn j 

where F := E + E^l"^ is the upper end of the window. Similarly, 

||i^(A + ii;)0||i.(,,) <C,2 / ||V(A + i?)^||2^C2^Qc, / (A + £;)0,(-A)(A + £;)0, 

Jn Jn 

= Cl J2 c]E, {E-E,f< ClEF . (2.7) 

3 

Using Cauchy-Schwarz, the sum of the last two terms in (j2.4p is then bounded by 
2Ca\/EF. For the first term on the right of (|2.4p . we use Einstein notation [A, D] = 
dii{ajdj-) — ajdiij. After several steps, using integration by parts and 4i\dn = 0, we 
get 

' q^[A,D](^ - 2 / (5,aj)(9,0)5j0 + / (9,,aj)05j0 (2.8) 
Jo Jn 

The constants := sup^^^ ll''^(x)l|2 where the matrix A e M"^" has entries diUj, 
and C" :~ sup^gj^^-^j^ „ |Aaj(x)|, exist and are finite. Then (|2.8p is bounded by 
2C'^F + C"i^^/2. Adding all bounds on terms in (|2^ we get 

Wdn^Uan) < ^{Ca + C:)F + C^'/f , (2.9) 

which is bounded by a constant times i? for i? > 1. □ 

Proof of Theorem M.SX Consider the coefficient vector c := {cj} g appearing 
in (|2.ip . where N is the number of eigenvalues (counting multiplicity) in the spectral 
window. Define the linear operator T : — > L^{dil) by 

Tc^^c.V, (2.10) 
j 

Lemma O states that \\T\\i2^L2^on) < (CnE)^/^. Thus ||Tr*||i2(ao) < CnE. But 
TT* is the operator in the statement of Theorem 11.31 which completes its proof. 



^The computation, involving a total of three derivatives, is justified for our class of domains, since 
Dirichlet eigenfunctions are in H'^^^{Q) for any Lipschitz O; see |13| . Theorem B, pl64. Rellich-type 
computations are also justified on Lipschitz domains in [1]. 
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3. Relating tension to a boundary operator. In this section, we sliow, fol- 
lowing Barnett that the tension t[u\ is related to the operator norm of a natural 
boundary operator. 

For E a non-eigenvalue of il, let IC{E) : L^(dil) L^(il) be the solution operator 
(Poisson kernel) for the interior Dirichlet boundary- value problem, 

(A + E)u = infl (3.1) 
u = f ondn , (3.2) 

that is, u = JCf. (For existence and uniqueness for data on a Lipschitz boundary 
see for example [23l Thm. 4.25].) Since the eigenbasis is complete in L^(n), we may 
write u = X^jli '^j'Pj- We evaluate each Cj by applying Green's 2nd identity, 

{E - Ej){d^j,u)L2(n) = / {u/\(j}j - (j)j^u) = [ {fij, - dp,dnu)ds , (3.3) 
Jn Jon 

thus Cj = {ipj, f) /{E — Ej). The solution operator may therefore be written as a sum 
of rank-1 operators, 

^(^)-f:|%^- (3-4) 

By the definition (|1.4I) we have, now for any u satisfying (A + E)u = in £7, that 
^[w] "'^ < ||/C(£')||. Since ||/C*/C|| = then by defining the boundary operator in 

A{E) := JC{E)* IC{E) , (3.5) 

we have an estimate on the tension that will be the main tool in our analysis, 

t[u]-^ < \\A{E)\\ . (3.6) 

Inserting p.4p into p.5p and using orthogonality (or see [H Sec. 3.1]), we have that 
A also may be written as the sum of rank-1 operators, 

This sum is conditionally convergent: the sum of the operator norm of each term 
diverges. For instance, for n = 2, Weyl's law [171 Ch. 11] states that the density of 
eigenvalues Ej is asymptotically constant, but since Hf/'jlP = ^(^j) the sum of norms 
is logarithmically divergent; for n > 2 the divergence is worse. Despite this, we have 
the following, which improves upon the results of [4]. 

Lemma 3.1. Let Q, CR", n>2, satisfy Condition[l\ and let E > 0. Then 

converges in the norm operator topology. Furthermore, the limit operator A{E) is 
compact in L^{drt). 

Proof. This follows immediately from (|4.8|) in the proof of Lemma 14.21 below, 
which shows that the tail of the sum in (j3.7p has vanishing operator norm. A is 
therefore also the norm limit of a sequence of finite-rank operators. □ 

7 



4. Proof of Theorems 11.11 and 11.21 In the previous section we related tension 
to the norm of a boundary operator which itself can be written as a sum involving 
mode boundary functions. Here we place upper bounds on in order to prove 

Theorems 11.11 and 11.21 Firstly we note that when E is an eigenvalue, Theorem 11.11 is 
trivially satisfied, since i[u,-„i„] — 0. When £' is a non-eigenvalue, formula (|3.7p enables 
us to split up contributions from different parts of the Dirichlet spectrum. 



A{E) = A„_(i;) + A^^{E) + A^4E) (4.1) 



where 



A^UE)= y -M^^ (4.2) 

El2<Ej<.1E,\Ej-E\>E^I'^ ^ ■" 

^-(^)" 2. (^E-E,y^ ^ {E-E,f ^^-^^ 

Ej<E/2 ^ ■" Ej>2E ^ ^' 

It is sufhcient (due to the operator triangle inequality) to bound the norms of these 
three terms independently. We first tackle the "far" and "tail" terms. 
Lemma 4.1. There is a constant C dependent only on such that 

\\A^AE)\\ < C forallE>l (4.5) 



Proof. For any E > 1, consider the spectral interval Im '■— [E + mE^/^,E + (m + 
1)E^/^]. For any such interval lying in [E/2, 2E] we may apply Theorem II. 3[ with E 
replaced by at most 2E, to bound || J2e ei V'j(V'ii)ll by 2CfiE. For terms in (|4.3p 
associated with this interval, the denominators are no less than m?E. Thus 

II V'j(V'jv) < ^ UQ) 

II ^ {E-E,Y - m2 ^ ' 

Covering [iJ+E'^/^, 2E] by summing over m = 1, 2, . . . gives a constant, since ^ m^^ = 
TT^/G. The same argument applies for intervals covering [E/2, E — E^/"^]. □ 
Lemma 4.2. There is a constant C dependent only on f2 such that 

\\A,^^x{E)\\ < CE-^I^ forallE>l (4.7) 



Proof. Consider a spectral interval := [2"^E, 2™~^^E]. We may cover this with 
at most 2™/2~ii?i/2 + 1 windows of half-width at most 2"^^^E^/^; for each of these 
windows Theorem 11.31 applies to bound || X^BjG/™ V'i(^j7)ll by Cn2"'+^i?. For each 
Ej e Im, the denominator is no smaller than (2™~^£')2. Thus 



, ^ (1^1 - (2'"/^-^i?^/^+l)P^=C^n(2-"/^-^i?-V2+2-"+i^-i) 

(4.8) 



The infinite sum over m — 1,2, . . . gives 



I ^ (E 

Ej>2E ^ 



< Cc. 



4(V2-1) 



+ 2E-'^ ] < CE-^I"^ for all E > \. (4.9) 



We treat the interval (0, £'/2) similarly, using a sequence of intervals J™ := [2 ^i?, 2 
Each such interval may be covered by at most 2"(™+'^)/^iJ^/^ + 1 windows of half- 
width 2-('"+i)/2£'i/2_ For each Ej e J„, the denominator is no smaller than E"^ j^. 
In a similar manner as before, the operator norm of the partial sum associated with 
J,„ is then 0{2-'"^E-^l'^), thus the infinite sum over m is 0(E-^l'^). Note that The- 
orem [L3] does not apply for i? < 1, but that there are 0(1) such Ej values and each 
contributes OiE^^). This proves the Lemma. □ 

Proof of Theorem ll.ll Examining the "near" term (j4.2[) . we use Theorem 11.31 on 
the sum of numerators, and get a bound by taking the minimum denominator. 



U_(i?)|| < 



CqE 



for all > 1 



diE,ay 

Using this and the above Lemmas to sum the terms in (j4.ip gives 



d{E, a) 



C 



for all > 1 



(4.10) 



(4.11) 



From Lemma IB.li an upper bound on the distance to the spectrum, we see that the 
second term is bounded by at most a constant times the first, so may be absorbed 
into it to give 



< 



CE 



d{E,aY 



(4.12) 



for all > 1 



Combining this with (13. 6|) proves (11.81) , hence also the second inequality in (|1.9p . The 
first inequality in (|1.9p simply follows from the fact that, since A is a sum of positive 
operators, 



t\u„ 



\A{E)\\ > 



[E-E, 



d{E,ay 



where Ej is the eigenvalue closest to E. Using the lower bound 
[TO] this becomes 

d{E,a) > c^/E]t[u^,,-,] . 



(4.13) 
> cEj from 

(4.14) 



With a change of constant, Ej may be replaced here by E to give the first inequality 
in (|1.9p . since Lemma FB . 1 1 insures that Ej is relatively close to E. (The lemma is not 
useful for E less than some constant and Ej < E, but then the ratio E/Ej is still 
bounded by a constant because Ej > Ei). 

Proof of Theorem \l.Sl We next prove the eigenfunction error bound (fTTOl . first 
considering E a non-eigenvalue. We denote the boundary data by U := u|asi- From 
orthogonality, then using the formula for the Ci coefficients below (13. 3p . we get. 



[E-E,f 



< 



IE 

E.=^E, 



(E - E,y 



(4.15) 



^E] 
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The operator in the last expression is identical to (I3.7P except with the iJj-eigenspace 
terms omitted. Therefore, its norm may be bounded in the same way as that of 
A{E), the only difference being that the d{E,a) introduced in (|4.10l) is replaced by 
vaiWE^^Ej \E — Ei\ = \E — Ek\- Thus the bound analogous to (|4.12l) is 

II V- ^J-(V'jv) II < CE foraUK>l 

and inserting this and ||?7||L2(af2) = ^[u]||u||2,2(f2) = t[u] into (|4.15l) gives (ll.lOp . Fi- 
nally, if E is an eigenvalue, i.e. E = Ej, the solution operator (13. 4p is undefined, since 
a solution to (|3.ip - p.2p exists if and only if / is orthogonal to the normal derivative 
functions in the _B-eigenspace. This can be seen by applying Green's 2nd identity 
to (p, any function in the _B-eigenspace, and u, giving (i9„0, U) — 0. However, the 
solution coefficients Ci for which Ei ^ E are uniquely defined by the same formula as 
before. Thus (|4.15p and the rest of the proof carries through. 

Remark 4.1. Theorem \1.2\ is sharp, as can be seen in the following way: if u is 
such that t[u\ is close to t[u,ni„] (say, less than 2i[u,„i„]j, then we have, by combining 
(0 and (fTTUl) . 



^^m^m < c^-§^ • (4.16) 



Apart from the value of the constant, one cannot expect to do better than this. For 
example, if E is midway between Ej and Ek, then the error \\u — Uj||L2(f2) cannot be 
expected to be better than l/\/2. 

5. Star-shaped planar domains. The purpose of this section is to say some- 
thing stronger than Theorem 1 1.11 in the special case of star-shaped domains in n = 2. 
We take weighted boundary functions 

^pfis) (x(s) • n(s))9„(/.j(s) , s e dQ (5.1) 
and our boundary inner product as 

{f,g)s:= [ i^{s)-nis))-'f{s)g{s)ds (5.2) 
Jan 



hence norm := (/, /), and ts[u] :— ||[/||s/||M||L2(r2). The significance of the 

weight (x • n) is twofold: it is strictly positive for strictly star-shaped domains, and 
also turns the inner product in (|1.13|) into s, enabling us to benefit from 

pairwise quasi-orthogonality. The Rellich theorem ||V'j''''||s = states that, with 
this special weight, there is no fluctuation in the L^-norms of the boundary functions. 
As shown in [4 , the function is[u,-„i„] vs E has slope in the neighborhood 

of Ej (this arises from dominance of a single term in (|5.5p below) . Hence these slopes 
are predictable independently of the particular form of each mode (/>j . This enables us 
to get the following eigenvalue inclusion result analogous to Theorem ll.il 

Theorem 5.1. Let f2 C K.^ be a strictly star-shaped bounded domain with piece- 
wise smooth boundary. Then there are constants ci, C2, C3 depending only on fl, 
such that the following holds. Let E > 1, and suppose u is a non-trivial solution of 
{A + E)u = m C°°{n), with C2ts[u]'^ < I. LetF:=E + VE. Then, 



(5.3) 



For the Helmholtz solution u,„i„ minimizing ts[u\ at the given E, 



2{E - c^E^/^)ts[u^,J\ < d{E,(j) (5.4) 

Remark 5.1. In the limit of high frequency i? 3> 1 and small tension tg [u] <^ 
E"^^^, the right-hand side of (|5.3p and the left hand side of (|5.4p are both ^/2E{1 + 
o(l))ts. This proves that both the power of E and the constant \/2 are sharp. 

Remark 5.2. Notice that this theorem is not applicable for all E since there may 
be large spectral gaps where C2ts[u]'^ < 1 cannot be satisfied. Due to the numerator, 
it becomes far from optimal when tg[u] is 0{E~^/'^) or larger. In these respects it is 
less general than Theorem ll.ll even though it gives better bounds in the small tension 
limit. 

The main tool used in the proof of Theorem l5.1l is the pairwise quasi-orthogonality 
result, Theorem II. 4[ together with the Cotlar-Stein lemma, which we state here for 
the special case of self-adjoint operators: 

Lemma 5.2 (Cotlar-Stein [T21|311[II1). Let{Tj}j^j be a countable set of bounded 
self-adjoint operators, J C N. Then 



|pT,|<mj.xj;^||T.T,|| . 

jeJ ieJ 

Proof of Theorem I5.il The weighted equivalent of p.7p is the operator 

which, by analogy with p.6|) . satisfies 

ts[u]-^ < WA^'Hms . (5.6) 

The lower bound follows by analogy with (|4?T3l) - (|4?T4|) . using = 2Ej, 

and Ej > E — c^E^/'^ from Lemma IRT] 

Using the same splitting into "near" , "far" , and "tail" parts as in Section IH we 
can bound the norm of the "near" part in a new way, as follows. 

Lemma 5.3. There is a constant ci > depending only on 57 such that 

2F V2ciF 
d{E, erf d{E, a) 

The first term in this bound will arise simply from the single term in the sum (j5.5p 
with Ej nearest to E. The second term requires more work, as we now show. 

Proof Let J = {j : \Ej - E\ < E^/^}. Using Tj = j^^%4i in Lemma [5j gives 



< „t 2 +17^ forallE>l 



(E-E, 

Applying quasi-orthogonality fTheorem ll.4p for the inner product, and |! II? '^Ej, 
and separating diagonal (i = j) from off-diagonal terms, we get. 



\j^,{E-E,Y ^fif^E^E,Y+'^^^'fif\E-E,\2^^ \E-E^ ^'''^ 



2E, ^„ e]/' ^e]"\K-E 
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Here S is as in Theorem 11.41 The first term is bounded by 2F/d{E,a) . Using 
\E.i - Ej \ < \Ei - E\ + \E - Ej \ bounds the second term by 

Sr^ j^i/i f- \E-Ei 



\E-E^\f^j ' V \E-E, 

< V2FS\J\ + -J—) < (5.9) 

jeJ \\E - Ej\ d{E,a)J d(E,a) 

Recall Weyl's law for the asymptotic density of eigenvalues, which states that, for 
n = 2 and vol = 1 , 



N{E) #{j : E, < E} ^ ^E + R{E) , (5.10) 



where the remainder is R{E) = 0{\ E) (|29|; for the case of piecewise-smooth bound- 
ary see [31, Eq. (0.3)]). Since the remainder is bounded for small E, there is a constant 
Cw such that \R[E)\ < C^^/E for all E > I. Thus | J|, the number of terms in the 
"near" window, is bounded by 



\J\<{^ + 2C^)Vf. 



Inserting this into (15.91) proves the Lemma, and we may take ci = 25'(l/47r-|-2Cw)- □ 
Completion of the proof of Theorem \5.1i The proofs of analogously weighted 
versions of Lemmas 14.11 and 14.21 are unchanged. So we may combine them with 
Lemma 15.31 and (j5.6p to get, for some constant C2, 

ts[ur' < \\A('\E)\U < ■:;;^ + ^k^ + C2 forall£;>l 

Multiplying through by d{E,a-)'^ we solve the quadratic inequality for d{E,a), 
^ . cif /x/2 + ^clFy2 + 2{ts[u]-^ 

ts [u\ ^ - C2 

Using the subadditivity of the square- root completes the proof of (|5.3p . 



6. Discussion of explicit constants. For the practical application of Theo- 
rems 11.11 and 11.21 it is important to have an explicit value for the constant C (from 
the discussion after (j4.15p we notice that C in the two theorems is the same.) 

We now compute an explicit value of this C that holds for all i? > 1. Examining 
(j2.9p we see that a choice of constant in Lemma 12.11 a-nd hence Theorem 11.31 that 
holds for all > 1 is Cn = 4(Ca + C'^) + V2C'^. To compute this we need sup norms 
of the value, and first and second derivative, of a vector field a as in Lemma [2?2] 

The proof of Lemma 12.21 shows such a construction; the values will depend on 
the size of the vectors a^,. and the choice of partition of unity used to cover fi. The 
vectors a^^ will be large (order 1/e) if has corners with angles less than e or greater 
than 27r — e. We note that a numerical procedure for this construction could be useful. 

In some special simpler prescription for the vector field can be given: 

• For strictly star-shaped domains in M", we may choose a = x/ infaf2(x • n), 
which gives Ca = supg^^ (x- n)/ infao(x- n), C"^ = 1/ infao(x- n), and C" = 0. 
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• For a domain with boundary, let (5 > be the largest number such that 
for each xq G dfl, a ball of radius S can be placed within 57 so as to be 
tangent to d^l at Xq. We may then choose a = (1 — r/S)'^nr, for r < (5, a = 
otherwise, where the coordinate r is the distance from dfl, and n,. is the 
unit vector in the local decreasing r direction. This gives constants Ca = 1 
and = 2/6. C'^ depends on 5 and an upper bound on the rate of change 
of surface curvature. (Also note that a slight modification of the proof of 
Theorem 11.31 would allow estimation purely in terms of Ca and , but with 
a doubling of the numerical constants). 
Summing the terms (j4.6p above and below E we have that the constant in 
Lemma [4.11 is 27r^Cn/3. Similarly, using (|4.9I) and its equivalent for (0,i?/2) gives 
the constant in Lemma l42l as [) + 6) < 7Cn- Summing these two 

constants gives a constant C in ()4.1ip as 14Cn. A choice of constant in (|4.12p is then 
Csi + 14Cn max[i5^, Cj], where from Appendix IB] we have Cd = and the max 

accounts for the case 1 < E < Ei. Finally, the constant in (|1.8p is the square- root of 
this, C = VCn(l + 14max[£:i,4]£;i). 

Requiring that the above estimates hold for all > 1 caused non-optimality in 
the choice of constant. It is more sensible in high frequency applications to use a 
better constant which is approached for ^ 1, and small tension t ^ 1. We now 
give this explicitly. In this limit, in (j2.9|) . F tends to E, and we drop lower-order 
terms to get Cn = 2{Ca + C^), which in the star-shaped case is 

Cn = 2i±£l^PMill_fl) for S > 1, 17 star-shaped . (6.1) 
mlaj2(x • n) 

If tension is small (i.e. E is not in a large spectral gap), the second term in (|4.12p is 
negligible, so we may approximate the constant in (jl.Sp as 

C = y/C^ for £; > 1, t < 1 . (6.2) 

Remark 6.1. The limiting constant (16.21) does not reflect the limiting slopes of 
the graph i[u,-„i„] vs E near eigenvalues. These slopes are known to he l/WipjW^, 
which is hounded by (2C^i?j)~^ l^d]- 

We end by discussing the constants ci and C2 in Theorem 15. II Constant C2 may 
be estimated easily, as above, using the weighted versions of Lemmas 14.11 and 14.21 In 
the proof of Lemma 15.31 ci involves the Weyl constant Cw; we know of no explicit 
estimates for Cw in the literature (the closest we know are estimates of the form 
|i?(£^)| < C^/ElnE with explicit constants [TU]). However, these constants are 
effectively irrelevant for practical purposes, when E ^ 1 and t <^ E~^^^, since in 
these limits, one may replace (|5.3p by d{E,a) < \j2Ets \u\ and still have an error 
bound very close to that given by the full expression. 

7. Numerical example. In Fig. II. lb we show a planar nonconvex domain given 
by the radial function r(6') = 1-1- 0.3cos[3(6' -f- 0.2 sin 0)]. The domain is star-shaped 
and smooth (we will not address numerical issues raised by corners here; see p!6l [T4l 
El I15| [31 [3.) For high-frequency eigenvalue problems, a convenient computational 
basis of Helmholtz solutions are 'method of fundamental solutions' basis functions 
^„(x) = yo(v^|x — yn|)7 where Yq is the irregular Bessel function of order zero, and 
{yn}^=i ^ of 'charge points' in M^\17. The latter were chosen by a displacement 
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of the boundary parametrization x(0), < 9 < 27r, in the imaginary direction (see 
[B]); specifically y„ = x(27rn/iV - 0.025i). 

We compute the data plotted in Fig. ll.lb . b as follows. At each E, i[u„i„] is given 
by the square-root of the smallest generalized eigenvalue of a generalized eigenvalue 
problem (GEVP) involving N x N symmetric real dense matrices F and G (the basis 
representations of the boundary and interior norms respectively.) Both matrices are 
evaluated using M-point periodic trapezoidal quadrature in 9, that is, quadrature 
points Xm — x(27rm/M), m = 1, . . . , M, and weights Wm = 27r|x'(27rm/Af)|/M. For 
instance, F ^ P*P, where P G M'*^^^ has elements 

Fmn — ^/^m (x^m ) : ('^•-^) 

and G is similarly found [4|, Sec. 4.1] using P and the matrices P*^^^ and P^^^ whose 
entries are the xi- and a;2-derivatives of those in P. Since the GEVP is numerically 
singular, regularization was first performed, similarly to [5J Sec. 6], by restricting to an 
orthonormal basis for the numerical column space of [P; P*^^); P^^^] given by the left 
singular vectors with singular values at least lO"^"* times the largest singular value. 

For low frequencies (Fig. II. ik ). 8-digit accuracy requires N = 100 basis functions 
and M — 200 quadrature points. For higher frequencies corresponding to 40 wave- 
lengths across the domain (Fig. II. lb . c), it requires TV = 400 and M — 500, and the 
above GEVP procedure takes 3 seconds per E valueH Very small (< 10^^) tensions 
cannot be found this way, and instead are best approximated via the GSVD [8]: the 
optimal tension at a given E is the lowest generalized singular value of the matrix 
pair (P, Q), where matrix Q has entries 

^ /x„i ■ Jim I I In n\ 

Qmn = d ^ -^(^™) ' ('-2) 

where is the normal at Xm, and regularization as before. Note that G ~ Q*Q well 
approximates the interior norm in the subspace with zero Dirichlet data, due to the 
Rellich formula (case i = j oi Theoreni ll.4[) . 

Any single- variable function minimization algorithm may then be used to search 
for a local minimum of t[M„,i„] vs E; we prefer iterated fitting of a parabola to t[M,„i„]^ 
at three nearby E values, which converges typically in 5 iterations. Using this with 
the GSVD (with N — 500, M = 700, i.e. 6 points per wavelength on dfl, and taking 
8 seconds per iteration), we find the tension 

= 2.2 X 10"^2 at E= 10005.0213579739 . (7.3) 

This is shown by the dot in Fig. II. lb . The GSVD right singular vector gives the basis 
coefficients of the corresponding trial function M„,i„, which is plotted in Fig. II. lb (this 
took 34 seconds to evaluate on a square grid of size 0.005, i.e. 1.3 x 10^ points.) 

Armed with datapoint (17. 3p . what can we deduce about Dirichlet eigenpairs of 
using our new theorems, and how much better are they than previous results? 

______ -j^ /2 

The constant in the Moler-Payne bound (jl.Sp is Cmp — Qi where qi is the lowest 
eigenvalue of a Stekloff eigenproblem on [501 (2-11)]. Since ft is star-shaped, the 

1 /2 

bound qi > E{ infan(x • n)/2supgf2(x • n) from [5D1 Table I] applies, giving Cmp — 
1.31 as a valid choice. Thus (|1.5p states that there is an eigenvalue Ej a distance no 



■^AU computation times are reported for a laptop with 2GHz Intel Core Duo processor and 2GB 
RAM, running MATLAB 2008a on a linux kernel. 
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more than 2.9 x 10^* from the above E. On the other hand, (|6.2p and (j6.1|) give the 
constant in Theorem 11.11 as C = 2.9. Applying the theorem gives a distance from 
the spectrum of no more than 6.3 x 10"^". Taking the smaU-tension hmit of the 
star-shaped planar result (15. 3p . and recomputing the weighted tension is['u,„i„] from 
Section [5l we get an even smaller distance of 3.5 x 10^^°, that is, an error of ±3 in 
the last digit of (|7.3|) . The latter is a 80- fold improvement over Moler-Payne. (Also 
see [4] for an example at higher frequency with 3 digits of improvement). 

How good an approximation is u„,i„ to the eigenfunction (jjj? Using the observation 
that the next nearest eigenvalue is Ek = 10007.339 • • • , the eigenfunction bound of 
Moler-Payne [24] gives an L^-error of 1.2 x 10~®. With the same data, using the 
constant C above. Theorem 1 1 . 21 gives an L^-error of 2.7 x 10^^", a 50-fold improvement 
over that achievable with previously known theorems. 

8. Conclusions. We have improved, by a factor of the wavenumber, the Moler- 
Payne bounds on Dirichlet eigenvalues and eigenfunctions which have been the stan- 
dard for the last 40 years. This makes rigorous the conjectures based on numeri- 
cal observations in [4]. We expect this to be useful since high-frequency wave and 
eigenvalue calculations are finding more applications in recent years. Of independent 
interest is a new quasi-orthogonality result in a spectral window (Theorem II. 3p . 

For numerical utility, throughout we have been explicit with constants, and have 
specified a lower bound on E for which the estimates hold (this being stronger than 
merely a 'big-0' asymptotic estimate). For star-shaped domains we strengthened the 
inclusion bounds (Theorem l5.ip . achieving a sharp power of E and sharp constant, in 
the limit of small tension, when tension is weighted by a special geometric function. 
This weight allowed pairwise quasi-orthogonality to be used, but since an upper bound 
for the number of eigenvalues in a ^/E window is needed, this is only useful for n = 2 
(planar domains). 

We applied our theorems to a numerical example, enabling close to 14 digits 
accuracy in a high-lying eigenvalue, and 10 digits in the eigenfunction. Both are two 
digits beyond what could be claimed with previously-known theorems. 

Our estimate C^/Et[u] on the distance to the spectrum is sharp (up to constants) 
if the tension t[u] (or ts[u]) is comparable to <:[u„i„] (fs[u„i„]). However, numerically, 
one generally has access to other properties of u (e.g. its normal derivative) which can 
give more detailed information about the spectrum. For example, the powerful 'scaling 
method' [36l [3] is able to locate many eigenvalues using an operator computed at a 
single energy E. In another direction. Still [33l Thm. 4] obtains improved inclusion 
bounds when the approximate eigenvalue is equal to the Rayleigh quotient; in this 
case, the bound is proportional to tlu]"^, but scaling as E^ for large energy. 

An open problem with practical benefits is the generalization of these results to 
Neumann and Robin boundary conditions, and to multiple subdomains with different 
trial functions on each subdomain and least-square errors on artificial boundaries 
[MKHIIT] (these are known as TrefFtz or non-polynomial finite element methods). 
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Appendix A. Proof of Lemma 12.21 

Proof. For every boundary point x, we can find a constant vector field a.^ liaving 
tfie property (|2.3|) in a neiglibourhood Ux of x. (Take a multiple of the vector field 
dz^ in the Euclidean coordinate system used in (jl.7l) .) By compactness we can find 
a finite number of such neighbourhoods Ui — Uxi, i — 1, . . . , N covering dil. We can 
add to this collection of open sets one additional set Uq, whose closure does not meet 
dfl, yielding an open cover of fl. Let i = . . . N he a smooth partition of unity 
subordinate to this open cover. Then 

N 

is a vector field with the required property. □ 

Appendix B. Upper bound on distance to spectrum. 

Lemma B.l. Let f2 C M" he a bounded domain. Let Ei he the lowest Dirichlet 
eigenvalue ofil. Then for any E > Ei, 

d{E,a) < CdE^/^ , Cd^ 2^1. 

Remark B.l. The result becomes interesting only for E > 2(1 + •\/2)£'i . Bounds 
on El exist as follows. Lf f2 contains a Euclidean ball of radius r, then Ei is less 
than or equal to Ei{B(0,r)), which is equal to i/'"^; where j„/2-i,i first 

positive zero of the Bessel function jn/2-i- Eor n = 2 we have jo,i = 2.4048... and 
for n — 3 we have ii/2,1 — i"- Also, Ei is greater than or equal to the first eigenvalue 
of the ball having the same n-volume as SI, by the Faber-Krahn inequality J27jj. 

Proof. Choose a wavevector k e M" with |kp = E — Ei, and consider the trial 
function m : f2 — >■ C defined by u(x) :— 0i(x)e*''''', where (pi is the normalized first 
Dirichlet eigenmode of SI with eigenvalue Ei. We calculate, 

(A + E)u = 2ik ■ V(f>ie'^-'' . 

Since u is in the domain of fl, has norm ||u||^2(q) = 1, and 

\\{A + E)u\\^ = 4 / |k- V0ipdx < 4{E- Ei)Ei < AEEi, 
Jn 

we see that u is an 0{E^/'^) quasiniode. On the other hand, writing u 
we have 

||(A + E)uf = II ^ a,{E - EM\^ = " 
3 3 

>d{E,afJ2\^3\'^d{E,af. 
3 

Combining (|B.ip and (jB.2|) completes the proof. □ 



(B.l) 
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